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ABSTRACT 



Monthly mean values or temperature output for a t.en-year 
period (1969-78) from a three-dimensional primitive equation 
ocean model are used to compute isotherm displacements at 
eight selected grid points in the North Pacific Ocean. The 
ocean model input parameters are climatological radiation, 
clouds, surface air tempejsture, surface humidity and six- 
hourly synoptic winds taken from Fleet Numerical 
Oceanography Center analyses. The first five baroclinic 
modes are computed from the mean temperature profiles at the 
eight selected model grid points and these modes are then 
fit in the time domain to the isctherra displacements. 

The fit shews evidence of baroclinic free Rossby waves 
in the time demain at the selected grid points south of 40N 
and indicates the predominance of the first baroclinic mode. 
At the selected grid points north of 40N, the results are 
inconclusive. This is due to the longer periods required at 
the northern latitudes, the neglect of salinity, and masking 
by a strong seasonal signal caused by surface heating and 
cooling at higher latitudes. 



4 



1 




TABLE 0? CONTENTS 



I. INTFCDUCIICN 10 

II. MODEL DESCRIPTION AND ANALYSIS PROCEDORE 12 

A. MODEL DESCRIPTION 12 

B. ANALYSIS PROCEDORE 14 

1. Governing Equations for Baroclinic Rcssby 

Waves 14 

2. Determination of N^ and The Isotherm 

Displacement Field 19 

3. The Solution of The Vertical Problem ... 28 

III. RESULTS 46 

IV. CONCIUSICNS AND RECCMMENDA TIC NS 65 

LIST OF REFERENCES 68 

INITIAL DISTHIBOTION LIST 69 



5 



LIST OF FIGURES 



Figure 1. Horizontal Domain of Ocean Model 13 

Figure 2. Profiles of (N^/^at Latitudes 18. 3N and 24. 4N 23 

Figure 3. Profiles of (N^)'^^at Latitudes 30. 5N and 36. 6N 24 

Figure 4. Profiles cf (lP)’'^^at Latitudes 42. 7 n and 48. 8N 25 

Figure 5. Profile of Latitude 54. 8N 26 

Figure 6. Isotherm Displacement at 30. 5N 139. 1M for 

Levels 1 to 6 29 

Figure 7. Isctherm Displacement at 30. 5N 139. 1W for 

levels 7 to 12 30 

Figure 8. Isctherm Displacement at 30. 5N 139. 1W for 

Levels 13 to 19 31 

Figure 9. Eiqencurves and Eigenvalues for Latitude 

30. 5N 34 

Figure 10. Ncrmalized Eigenfunctions at 18. 3N 37 

Figure 11. Ncrmalized Eigenfunctions at 24. 4N 38 

Figure 12. Normalized Eigenfunctions at 30. 5N 39 

Figure 13. Normalized Eigenfunctions at 36. 6N 40 

Figure 14. Ncrmalized Eigenfunctions at 42. 7N 41 

Figure 15. Ncrmalized Eigenfunctions at 48. 8N 42 

Figure 16. Ncrmalized Eigenfunctions an 54. 8N 43 

Figure 17. (t) and DEV(t) at 18. 3N 155. 9M 47 

Figure 18. (t) and DEV (t) at 24. 4N 155. 9W 48 

Figure 19. (t) and DEV (t) at 30. 5N 155. 9W 49 

Figure 20. A„ (t) and DEV (t) at 36. 6N 155. 9W 50 



6 



51 



Figure 21. A„ (t) and DE7(t) at 42. 7N 155. 9W 

Figure 22. A„ (t) and DEV (t) at 48. 6N 155. 9« 52 

Figure 23. A„ (t) and DEV (t) at 54. 6N 155. 9W 53 

Figure 24. A„ (t) and DEV(t) at 30. 5N 139. 1W 54 

Figure 25. A„ (t) , and DEV (t) for wsarhsr station 

November 56 

Figure 26. Cutoff Period versus Latitude 58 

Figure 27. Locus of Wave Number Vectors for Mode 1 at 

18. 3N 61 

Figure 28. Locus of Wave Number Vectors for Mode 1 at 

36. 6N 62 



7 



LIST OF TABLES 



TABLE I. 
TABLE II. 



Reciprocal Eigen val ues (m/sac) for Modes 1 to 



Mean Deviaticn {%) at Selected Grid Poinrs . 57 



8 



I 




'm 



ACKNCWLEDG EMENT 



The author expresses his sincere gratitude zc Or. Roberz 
Haney for net only his contributions as thesis advisor but 
also for providing the model simulated temperature data, for 
helping define the study and for providing recommendations, 
both editorial and theoretical. Without his help the devel- 
opment of this thesis would have been impossible. 

The author additionally expresses his gratitude to Dr. 
Andrew Willmott fer providing not only editorial recommenda- 
tions but also guidance pertaining to the mathematical and 
theoretical develcpment of the baroclinic Rossby wave equa- 
tion. Appreciation is also expressed to Dr. Rudolf 

Pr eisendcr f er for his guidance regarding the numerical solu- 
tion of the governing equation for baroclinic Rossby waves. 
Finally, the author expresses his deepest gratitude to his 
wife, Sheryl, and son, Michael, for their support and 
pa tience . 

The author utilized the computing facilities of the 
W. B. Church Computer Center at the Naval Postgraduate 

Schocl while ccnducting this study. 



9 



I. 



INTROD UCTION 



Oceanic general circulation models have been under 
development for the past decade and a half. Evaluation of 
these newly developed models was initially hindered by a 
lack of adequate observational data from ocean regions. 
However, ever the past eight years high quality oceanic data 
and data analysis studies from portions of the ocean have 
been made available to the large-scale ocean modeller. Of 
particular imccrtance and interest in this study has been 
the analyses cf temperature data resulting from the "Morth 
Pacific Experinent” (NORPAX) begun during the last decade. 

Eased upon an analysis of TRANSPAC data (a sub-program 
of NCEPAX confined to a region 30-50N and 170E - 150W), Kang 
and Magaard (1980) deteriined that annual temperature fluc- 
tuations in the area 30-40N, 160E - 130W consisted of first 

mods Rossby waves with random phase traveling northwest with 
a wavelength cf 300 km. For other latitudes of the Jforth 
Pacific Ocean, Price and Maqaard (1980) determined that 
first mode baroclinic Rossby waves with dominant period of 
seven to ten years existed in the region 40-50N. It was 
also determined that Rossby waves of annual period did not 
play a role in the latitude belt 40-50N in the North Pacific 
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Ocean . 



Magaard (1983) # ir. a paper discussing baroclinic 
Rossty wav€ energetics, summarizes these earlier studies by 
srating that first mode baroclinic Rossby waves dominate the 
internal temperarure f lucruatio ns in the region 20-30N in 
the North Pacific Ocean fcr a whole range of periods between 
five months and ren years. Between 30-40N only annual or 
near-annual first mode baroclinic Rossby waves are evident. 
It is the results of these studies that have motivated this 
investigation, which is esssntually an evaluation of whether 
a general circulation ocean model, driven by observed synop- 
tic winds, can simulate the above characteristics of the 
North Pacific Ocean Rossby wave field. Following a proce- 
dure outlined by Emery and Magaard (1976), this study exam- 
ines the possible existence of baroclinic Rossby waves 
inferred frcm temperature fields produced by the ocean 
model . 



1 1 



II 



• ^SPEL DESCRIPT ION A NAL YSIS PROCEDURE 

A. HCDEL DESCRIPTION 

The mcdel used in this study is a three-dimensional 
primitive equation model developed at the Naval Postgraduate 
School and described by Haney (1980) . A flat bottom and 
rigid lid are prescribed. The model is driven by seasonally 
varying (monthly climatology) solar radiation, clouds, sur- 
face air temperature, and surface relative humidity. The 
surface winds utilized are taken from six-hourly synoptic 
analyses, 1969-1978, made at the Fleet Numerical 
Oceanography Center, Monterey, California. Salinity is neg- 
lected. The model domain shown in Figure 1 is a rectangular 
basin extending from the Equator to 65N and from 145E 
125W. The horizontal grid spacing is approximately two 
degrees in latitude and three degrees in longitude. The 
model is 4000 m deep with 20 irregularly spaced levels, 14 
of which are above 1000 m. The boundary conditions pre- 
scribed at the bottom and sides are those of an enclosed 
thermally insulated basin. 
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Horizontal doiain of temperature field used for 
shown by enclosed region. Triangles mark 
selected for Rcssby wave analysis. 



analysis 
grid poi: 



Figure 1. Horizontal Domain of Ocean Model 



The model also includes a dynamic adjustment mechanism 
as described by Adamec, et. aj.. (1981). Essentially, the 
mechanism assures that if the local gradient Richardson num- 
ber, Hi, falls below a critical value (0.25) then the temp- 
erature and the velocity components at the two adjacent 
levels are mixed vertically so that heat and momentum are 
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cons€rv6d, tte aiixing ratios for hear and momentum are 
equal, and after mixing Ri equals the critical value. 

The model was run for a 30-year period, with the first 
two decades considered as a ”spin-up" period and the third 
decade, 1969-78, considered as a simulation of actual condi- 
tions. The purpcse of this thesis is to analyse the temper- 
ature variability simulated by the model during this last 
decade. The ccal is to see if the variability is consistent 
with the dispersion relation for baroclinic Rossby waves. 

B. ANALYSIS ESOCEOaRE 

"I ♦ Gove rni ng Eg uati cns for B aroclini c Ros sby Waves 

The governing equations for baroclinic Rossby waves 
are derived from the frictionless momentum equations, the 
hydrostatic balance, and the continuity equation for an 
incompressible fluid in the following form: 

+ (V*V)u - ^ 

V.J. + (V*V)v + fu + i Py = 0 

V*V = 0 

+ pg = 0 

+ (V*V)p = 0 



(la) 

(lb) 

(l c) 

(l d) 
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where V is the velocity vector, u, v and w are the velocity 
compcnents in the eastward (x) , northward (y) , and upward 
(z) direction, respectively, p is the density, p is the 
pressure, g is the acceleration due to gravity, f is the 

Coriolis parameter, t is time and V= *§t) • sub- 

scripts denote partial differentiation with respect to the 
subscripted variable. 

The vcrticity equation is formed using equations 
(la) and (1b). Subsequent application of the quasi- 
geostrophic approximation and linearization leads to the 
following fcrm of the vorticity equation; 



f ^ i 

^XX 



I ^ 





( 2 ) 



where p' and w' are the perturbation pressure and vertical 
velocity ccmpcnent, /8 = df/dy, p, is a constant reference 
density, and the other variables are as defined before. 
Additional assumptions used in the above linearization pro- 
cess are that the lateral shear of the basic flow is neg- 
lected and the basic state horizontal velocity is at rest 
(i.e., u = V = 0) . Prom this assumption, equation (1c) and 
the boundary condition that w = 0 at z = 0 (sea surface) it 
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follows that w = 0. Equation (2) contains two unknown vari- 
ables, p' and w*, therefore, an additional equation is 

required. 

Substituting from (Id) into (1e) after multiplying 
t>y “9/Po subsequent linearization yields the second 

equation , 

p * + p N^w' =0 (3a) 

^zt o 



or 

w' 



N^p 



■> 



(3b) 



where = ~(9/pjPz is ‘^-be 
assumed to be a function of 
using (3b) yields. 



(p' + p' ), +. 3p' + 

^xx ^yy t ^x ‘■j^2 



Brun t-Vaisala frequency squared, 
z only. Eliminating w' in (2) 



P'zth = “ 



(4) 



Using the separation of variables technique the perturbation 
pressure is written. 



p'(x,y, 2 ,t) = ’ 



(5) 
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where n is the mcde namber. Substituting (5) into equaricns 
(4) and (3fc) gives. 



_1 1 n 

P dz 

n N 



and 



] 



'' 'Pn>t " 6(P„) 



X 






= X , 
n ’ 



1 dP„ 



N 



2 dz 



PQ 

<Pn>t 



( 6 ) 



(7) 



2 2 

where V = (^2 ^^ 2 ) * -s a separation constant with 

dimensions of inverse velocity. Differentiating equation 
(6) with respect to z yields, 



d_ fi aSn + 1 2 = 0 

dz^ ‘,.2 dzl * dz 0 • 



( 8 ) 



Substituting- from equation (7) into equation (8) yields the 
governing equation or vertical equation. 



+ X ^ N^(z)w^ = 0 

dz^ 



(9a) 



subject to boundary conditions. 



w' = 0 at z = 0 (Sea Surface) 

w’ = 0 at z = -H (Sea Bottom) 



(9b) 
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Equation (6) written in rh€ form. 



+ /3(Pn)x = 0 



v2 ,2 



(10a) 



represents the Bcssby wave equation for which solutions of 
the form 



= EXP[i(kx <■ iy - ut)] , 



(lOb) 



are scuqht where k and Si are wavenumbers in the x and y 
direction, respectively and <r is the frequency. Substitut- 
ing from (10b) into (10a) yields a form of the dispersion 
relationship fcr baroclinic Rcssby waves, viz. , 



n 



- Sk - g(k^ + 
of^ 



( 11 ) 



Equation set (9) forms a Sturm-Liouville eigenvalue 
problem for the vertical motion, where A,, are the eigen- 
values and w* the corresponding eigenf unctions . A solution 
is sought for a finite number of modes, m. These modes cor- 
respond to the first m baroclinic Rossby wave modes. 

If N ( 2 ) were independent of depth at a fixed loca- 
tion and time, the solution for w* could be readily 
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determined analytically as rr igonomerric funcrcions of 

(n7TZ)/H (see for example Eedlosky (1979, pp. 356-360)), pro- 
vided 



n 



njT 

lai 



n = 1,2,3,... 



( 12 ) 



In cerzain instances an analytical solution to equation set 
(9) can be found provided N^z) is, cr approximates, a well 
defined function (see for example Willmott and Mysak(1980) 
who solve the vertical problem with N^z) as an exponential 
function cf z) . Realistic distributions of N^are depen- 
dent upon depth and generally not well defined functions of 
z, therefore the solution of equation set (9) is obtained 
numerically. 

2. Determin ati on of and The Isotherm Displacement 
Fi^d “ 

The Brunt-Vaisala frequency squared is, as defined 
earlier, equal to -(9/Po)Pz • Since the model neglects 

salinity, the simple state relation. 



p = p^Il - a(^ - T^)] , (13) 

is used to convert the Brunt-Vaisala frequency squared from 
a function cf the vertical mean density gradient into a 
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function of the vertical mean temperature gradient, where 
To is a ccnstant reference temperature and a is the coeffi- 
cient of thermal expansion. Differentiating equation (13) 
with respect tc z and substituting into the definition of 
gives. 



= ag . 



( 14 ) 



Eguaticn (14) is used to calculate from the mean 
monthly temperature fields output by the ocean model for the 
last ten years, 1969-78, cf the model run. Only the temper- 
ature field in the region 18.3 - 56. 9N and 159. IE - 125W 
shown as the enclosed box in Figure 1 is used in the calcu- 
lation of After completion of the numerical diff eren-^ ia- 

tion (right side of equation (14)), a ten-year' time mean, M 
at each grid pcint and at 19 vertical levels (mid-levels of 
the 20 model levels) ' is calculated. Finally N^is averaged 
over 28 east-west grid points, giving a versus depth pro- 
file for each latitude (20 north-south grid points) . Justi- 
fication fcr the east-west averaging is based upon little 
variability of with longitude over the region at a partic- 
ular vertical level. The variability of with latitude, on 
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the ether hand, is such thar north-south averaging is not 
justified. Price and Magaard ( 1980) also found that mean 
profiles varied much more with latitude than with longitude, 
and thus the numerical mcdel data are consistent with the 
corresponding data taken from the North Pacific Ocean. Pro- 
files of (N^ , averaged over longitude, are shown in 

Figures 2 to 5 for the latitudes of the grid points selected 
in this investigation (triangular symbols on Figure 1). The 
solid curves on each figure are drawn to the values of (N ^)’/2 
at 19 vertical levels determined by equation (14). The 
dashed curves are drawn to values of calculated at 

every five meters from z = 0 to z = -4000 m using interpola- 
tion. The interpolation is done numerically by means of 
Aitken's scheme of Lagrange interpolation described by 
Hildebrand (1956, pp. 49-50) using the values of the 

19 vertical levels and an assigned value of (N^/^ at -4000 m 
as tabular input. At latitude 48. 8N and 54. 8N the assigned 
bottom value for (N^ )’/^ was given the identical value calcu- 
lated at level 19 (2830 m) . At the other selected latitudes 
the assigned hettom value was determined by linear extrapo- 
lation using the levels below 10 00 m (levels 15 to 19) . 



21 



Th€ profiles shew a pronounced variability with 
latitude. Below approximately 900 m, (N^ is small (less 
than 0.001 1/sec) and nearly constant (model simulated 
isothermal condition in the deep ocean). At latitudes north 
of 47N this constancy is evident at depths below 150 m. The 
maximum value of (N^ )^^^^ (maximum vertical temperature gradi- 
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The second field to be computed using the model sim- 
ulated mean monthly temperature fields is the isotherm dis- 
placement. The isotherm displacement, D(z,t) , is the 
displacement of a monthly mean isotherm from its long-term 
average level. It is defined mathematically as -T'/('T )2 
where T’ is the temperature fluctuation, T’ = T - T. f and 
(T)j are the long-term (ten-year) averaged temperature and 
vertical temperature gradient, respectively. In the numeri- 
cal calculation of D(z,t) the assumption is made that the 
temperature field is linear in z between model levels. The 
isotherm displacement versus time curves at 19 vertical 
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Figure 2. Profiles of Latitudes 18. 3N and 2 <a.4N 
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Figure 4. Profiles of Latitudes 42. 7N ar.d 43. 3N 
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Figure 5. 
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Prefile of (N^ )’/2at Latitude 54. 8N 
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levels from the ten-year model simulation are shown in 
Figures 6,7, and 8 for the grid point at 30. 5N 139. 1W. 

By virtue of its definition, D(z,t) positive, in 
general, implies a temperature colder than the long-term 
mean. D(z,t) negative implies a temperature warmer than the 
long-term mean. The isotherm displacement curves, as 
depicted in Figures 6 to 8, reveal a pronounced seasonal 
signal (annual period) at levels above 71.5 m. At 100 m 
(level 6) and below the seasonal signal becomes less well 
defined and losses definition with depth. The seasonal sig- 
nal at the upper levels is primarily the model simulated 
mean monthly temperature fields response to the seasonally 
varying solar radiation provided as model input. Below 100 
m (Figures 7 and 8) a low freguency (long period) signal is 
evident. The isotherm displacement curves at the other 
selected grid points (not shown) , in general, show similar 
characteristics of an upper level annual cycle and low fre- 
quency signal at lower levels. The significant features, 
proceeding from south to north at longitude 155. 9W, are 
that, 1) the depth at which the seasonal signal remains 
distinguishable increases reaching 1760 m (level 17) at 
54. 8N, 2) the amplitude of the isotherm displacement at 
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mid-l'svels increases from a maximum of approximately 50 m ar 
18. 3N to 5000 m at 54. 8N. This increase in amplitude is a 

result of the near iscthemal conditions evident below 150 m 
at the northern latitudes (i.e. (T)^ is small thereby 

increasing the value of D( 2 ,t)). 

-• 2*2 € Solution cf The V ert ical P rob lem 

Equation set (9) is solved numerically using an 
Interactive Ordinary Differential Equation (lODE) solver 
which uses a fcurth-order Runge-Kutta integration technique. 
The eigenvalues, , are found by means of the '‘c-rcot 

methcd" described by Pr eisendorf er and Gonzalez ( 1 973) . In 
applying the "c-rcot method" equation (9a) is separated into 
two first-order differential equations, 

^ - 6 (15a) 

dz ^ 

^ N^(z)w‘ .• (15b) 

dz " " 

A trial value is chocsen for X, and equations (15a) and 
(15b) are integrated using lODE starting at the sea bottom, 
z = -H (4000 m) , with boundary values of w„' = 0 and 0= 1. 
As the integration proceeds the zeros of the function 
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ISOTHERM DISPLACEMENT 




1969 1970 1971 1972 1973 1974 1975 1976 1977 1978 



Level 1 = 5.0 m Level 2 = 16.5 m Level 3 = 33.0 ra 

Level 4 = 49.0 m Level 5 = 71.5 m Level 6 = 100.0 m 



Figure 6. I^ctherm Displacement ar 30. 5N 139. 1W for Level 
1 to 6 
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ISOTHERM DISPLACEMENT 
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l6V€l 19 is shown as dashed curve. 

Level 13 = 669.0 m Level 14 = 855.0 m Level 15 = 1089.5 m 
Level 16 = 1385.0 m Level 17 = 1760.0 m Level 18 = 2235.0 m 

Level 19 = 2830.0 m 



Figure 8. Isotherm Displacement 
13 to 19 



at 30. 5N 139. 1W for Levels 
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function c (z) - w»(z) - 0 (z)Z(z') ara noted. In the 

above expression Z(z') equals w*(z*)/ 0 (z*) where z* is the 
value of z az the bottom boundary which changes as zhe inte- 
gration continues. Since the boundary condition states that 
w' (-H) = w* (z') = Of then 2 (z ) must equal zero and the c 

function reduces to c (z) = w*(z). The eigenvalue problem 

thus reduces tc finding values of z where w*= 0 . 

Fcr example, for the p-ofile at 30. 5N and a 

trial value cf X= 2.5 sec/m, c(z) = w' (z) = 0 when z = 

-2U511I1, - 1 30 6 ni, -607m, -326m, and -112m. Varying \ by 0.1 
sec/m, noting when w'= 0 , and tracing out the points results 
in five curves (eigencur ves) as shown in Figure 9. The 
first five eigenvalues which satisfy the boundary condi- 
tions, equations ( 9 b), are found by the intersection cf the 
eigencurves with the z = 0 line. The first five eigenvalues 
for this example are shown to the right of the ordinate in 
Figure 9; namely, = 0.4735 sec/m, 0.8934 sec/m, 

X 3 = 1.3380 sec/m, X 4 = 1.7814 sec/m, Xg = 2. 1823 sec/m. 

Justification fcr limiting the calculation cf eigenvalues to 
the first five modes will be made apparent in the next chap- 
ter. An advantage in finding the eigenvalues in this manner 
is that the depth to the sea bottom need not be fixed. If 
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EI6ENCURVES FOR 30. 5N 



O 




Figure 9, Eigencurves and Eigenvalues for Latitude 30. 5N 

Having determined the first five eigenvalues by zh 
”c-rcot method”, the next step is to calculate the ccr 
respcnding eigenfunctions, , using lODE, Th 
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TABLE I 



Reel 


prccal Eig 


snvalues (m/s 


ec) for 


Modes 1 to 5 




LATITUDE 




MODE 


NUMBER 






1 


2 


3 - 


4 


5 


1 8 . 3N 


2.2578 


1.2442 


0.8505 


0.6274 


0. 5049 


24 .4N 


2 . 23 3 1 


1.1966 


0.8202 


0.6099 


0.4922 


30. 5N 


2. 11 19 


1 . 1193 


0.7474 


0.5614 


0.4582 


36 .6N 


1 .79 37 


C.9774 


0.6245 


0. 471 9 


0.3883 


42 .7N 


1.3289 


C.75 10 


0.4736 


0.3620 


0. 3077 


43 . 8 N 


0 . 92 9 9 


0.4831 


0. 3294 


0.2629 


0.2264 


54. 8 N 


0 . 55 1 1 


0.2598 


0.2153 


0. 1722 


0. 1356 



si genfunct ic 3 S are then rcrmal ized, following Emery and 
Magaard ( 1 976 ) , as 



Y„(z) = 

n o 



Wn'(z) 



/ ” 






-H 



( 16 ) 



where n represents tha mode number. Since equation ser (9) 
represents a Sturm-Licuville system, rhe eigenfunctions are 
orthogonal with respect to the weight function (z) over 
the vertical cclunm z = -H to z =0, viz. . 



O 

j N^(z)w^' (z)w^* (z)dz = 0 ; 

-H 



n 5^ m 



( 17 ) 
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Tharefore, ncrmalization of the eigenf unc~ions rasal-cs in 
the crthonornial sst,|y^| , viz. , 



where 6nm Kronscker delta function. A more detailed 

discussion ccnceining the orthogonality condition for tha 
eigenfunctions is given by Millmott and Mysak (1980) . The 
normalized eigenfunctions for the selected latitudes are 
shown in Figures 10 through 16. The characteristics of the 
eigenfunctions are related to the profile of (z) . At 
latitudes south of 40N (Figures 10 to 13) the nodal points 
of the normalized eigenfunctions are concentrated near the 
upper half of the column in response to the model simulated 
upper ocean maximum (maximum temperature gradient). At 
the northern latitudes (Figures 15 and 16), the normalized 
eigenfunctions are very nearly sinusoidal, particularly at 
the higher modes. This corresponds to the near constant 
profile of (N^ (Figures 4 and 5) below 150 m. 

The analysis procedure is concluded by fitting the 
normalized eigenfunctions to the simulated isotherm dis- 
placement, D(z,t). This is dons in the time domain by 



O 




( 18 ) 
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NORMRLIZED VERTICRL VELOCITY! 5EC/M^><l/2 ) 



- 20.0 - 10.0 0.0 10.0 20.0 




Mode 



Solid; Mode 
Mode 4 



2: Dashed; Mode 3: Long dash/Short 
Do-ted; Mode 5; Long dash/Dot 



Figure 10. 



Normalized 2 i genf unc tions at 13. 3N 



dash 
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NORMRLIZED VERTICRL VELOCITY( SEC/MhxI/2 ) 



- 20.0 - 10.0 0.0 10.0 20.0 




Mode 1: Solid: Mode 2: Dashed 
node 4 ; Dotted; 



Mode 3; Long dash/Short dash 
Mode 5; Long dash/Dot 



Figure 11. Ncrmalized Eigenfunctions at 24, 4N 
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NORMALIZED VERTICAL VEL0CITY(SEC/Mkk1/2] 



- 20.0 - 10.0 0.0 10.0 20.0 




Mode 



Sclid; Mcde 2: Dashed; Mode 3 
Mode 4: Dotted; Mode 5: 



Long dash/Short 
Long dash/Dot 



dash ; 



Figure 12. Normalized Eigenfunctions ar 30. 5N 
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NORMALIZED VERTICAL VEL0CITY(SEC/Mxxl/2 ) 



-30.0 -2CT.0 -10.0 0.0 10.0 20.0 30.0 




Mode 1: 



Solid: Mode 2: Dashed; Mode 3: Long dash/Short 
Mode 4: Dotted; Mode 5: Long dash/Dot 



Figure 13. Normalized Eigenfunctions at 36. 6N 



dash ; 
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NORMRLIZED VERTICRL VELOCITY ( 5EC/Mx><l/2 ) 



-30.0 -20.0 -10.0 0.0 10.0 20.0 30.0 




Mode 1: 
Figure 



Solid; Mode 2 
Mode 4 : 



14. Normalize 



: Dashed; Mode 3 
Dotted; Mode 5: 

d Eigenfunctions 



Long dash/Short 
Long dash/Dot 

an 42. 7N 



dash ; 



41 



NORNRLIZED VERTICfiL VEL0 CITY(SEC/Mhk1/2) 




Hode 



Solid; Mode 2; Dashed 
Mode >i : Dotted; 



Hode 3 
Mode 5: 



Long dash/Short 
Long dash/Dot 



dash; 



Figure 15. 



Ncrcnalized Eigenfunctions at 48. 8H 



42 



NORMRLIZED VERTICRL VEL0CITY(SEC/M><xl/2] 



-75.0 -50.0 -25.0 0.0 25.0 50.0 75.0 




Mode 1: 



Solid; Mode 2: Dashed; Mode 3 
Mode 4: Dotted; Mode 5: 



Long dash/Short 
Long dash/Dot 



dash 



Figure 16. Normalized Eigenfunctions at 54. 8N 
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making use cf the ccmpleteriess of the eigenfunctions to 
equate D(z,t) to the sum of the normalized eigenfunctions 
multiplied by an amplitude coefficient A^, ikl- > 

oo 

D(z,t) = ^ ^ • (19) 

n=l 



The coefficients, A„ (t) , are readily obtained by making use 
of the orthcncrmality condition. Multiplying each side of 
equation (19) by N^y^and then integrating from z = -H to z = 
0 results in 

O a. ® ' 

F(z)y (z)y (z)dzl . (20) 

n in 






(z)Yjj^( z)D( z ,t )dz 









-H 



n= 1 



/ 



-H 



Application cf the othogorality condition (18) to equation 
(20) yields 

O 

A (t) = / N^(z)y (z)D(z,t)dz . (21) 

n I n 

-H 

Equation (21) is used to compute A„(t) for the first 
five modes for each selected grid point. The devia'^non iron 
the modal fit, as defined by Emery and Sagaard (1976) ,eqn.(22). 
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is also caiculatsd. The deviati cr. term quan^iraxi vely shows 
how well the model simulated isotherm displacement flucxua- 
tions are fit by the first five baroclinic modes. A value 
of DEV(t) = 0 implies that the isoxherm displacemenx fluctu- 
ations can he fully explained by the first five baroclinic 
modes. 



o 



5 




# 



DEVCt) 



( 22 ) 



III. 



RESULTS 



The coefficients, A„(t), and the deviation, DEV(t) , cal- 
culated at each selected grid point are shown in Figures 17 
to 24. At lower latitudes (Figures 17 to 19 and Figure 24) 
the anplitude cf coefficient A, is generally larger than the 
coefficients A 2 tc Ag indicating the importance of the first 
mode at the southern latitudes. At latitude 24. 4M (Figure 
18) a seasonal signal is evident for coefficients A^ and Ag . 
The seasonal signal of the higher modes becomes more appar- 
ent with increasing latitude indicating that the seasonally 
varying model input parameters appear to be reflected in the 
higher modes at higher latitudes. At the northern most grid 
point, 54. 8N (Figure 23) the overall amplitude increases an 
order of magnitude. Additionally, the pattern of each coef- 
ficient appears nearly identical with only slight differ- 
ences in magnitude. 

The deviation term at latitudes 18. 3N and 24. 4N remains 
generally belcw 50% during most of the ten-year model run, 
although scae time periods are higher. North of 24. 4N there 
is a general increase over the ten-year period in the 
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LRT. 18.3N LONG. 155.9W 




Figure 17. (t) and DEV (t) at 18. 3N 155. 9W 
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LRT. 24. 4N LONG. 155. 9W 




Figure 18. (t) and DEV (t) at 24. 4N 155. 9M 
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Figure 19 



An (t) and DEV (t) ar 30. 5N 155. 9M 



LflT. 36. 6N LONG. 155. 9W 




Figure 20. (t) and DSV (t) at 36. 6N 155. 9W 
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gure 21. (t) and DEV (t) at 42, 7N 155. 9W 
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LRT. 48. 8N LONG. 155. 9W 




Figure 22. A„{t) and DEV (t) at 48. 8N 155. 9W 
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PERCENT (lliMK3/2)/S) ((Mx><3/2)/S) ((Mxx3/2)/S) ( (Mk>.-3/2)/S') I (Mkx3/2)/S1 

25.0 50.0 75.0 - 75.0 25.0 - 75.0 25.0 - 75.0 25.0 - 75.0 25.0 - 75.0 25.0 



LflT. 5-1. 8N LONG. 155. 9W 




Figure 23. (t) and DEV (t) at 54. 8N 155. 9W 
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LflT. 30. 5N LONG. 139. IW 




Figure 24. A„ <t) and DEV (t) az 30. 5S 139. IW 
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deviation curve with rhe exception of the two northern mcst 
grid points (Figures 22 and 23). As a comparison, the first 
five coefficients and the deviation term computed by Zmery 
and Bagaard (1976) using 48 months of data for weather sta- 
tion November are shown in Figure 25. The results shown in 
Figures 24 (model grid point closest to weather station 
November) and 25 compare reasonably well. Both show the 
largest amplitude of coefficeint A^ . The 48-racnth mean 
deviation computed for weather station November was 32% 
while the ten-year mean deviation is 38.2% for model grid 
point 30. 5N 139. 1H. 

To examine the relative contribution made by each baro- 
clinic mode tc fitting the simulated isotherm displacement 
fluctuations (egua tion ( 1 9) ) , the ten-year mean deviation 
(equation (22)) was calculated for each selected grid point 
adding one mode at a time. The results are given in Table 
II. Column 1 (under modes) gives the ten- year mean 
deviation when only the first baroclinic mode is fit to the 
isotherm displacement fluctuations. Column 2 gives the 
ten-year mean deviation when modes 1 and 2 are fitted to the 
isotherm displacement fluctuations and so on. The tabular 
values show the importance of mode 1 at the selected grid 
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1366 



1967 



1968 



1969 




Figure 25. An (t) , and DEV(t) for weather station November 
(from Emery and Hagaard, 1976) 

points south cf 30. 5N (i.e., over 40% of the deviation is 
accounted for by rhe firs- baroclinic mode alone) . The ccn- 
tribution made by the fifth barcclinic mode in reducing the 
mean deviation is less than 6 % at the three southern most 
grid points, and, except for the northern mosr grid point, 
rs less than elsewhere. Based upon these results, m rs 

apparent that adding higher baroclinic modes (modes greater 
than 5) wculd add little tc reducing the deviation term com- 
pared to mode number one. 
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lABLE II 





Mean 


Deviation 


(^) at 


Selected Gr 


id Points 




GRID POINT 






MOD^ 










1 


1+2 


1 + 2+3 


1+2+3+4 


1 +2+3+4+5 


1 8. 3N 


155. 9W 


56.3 


46. 0 


36. 4 


28. 6 


26.8 


24. 4N 


155. 9W 


54. 1 


45. 9 


39.0 


35. 0 


31.6 


30. 5N 


155. 9W 


57.4 


54. 2 


50. 9 


46.9 


4 1.0 


36. 6N 


155. 9W 


65.2 


55.7 


53. 3 


45.3 


38.8 


4 2. 7N 


155. 9S 


96 . 2 


88. 1 


60.7 


54. 2 


46. 2 


48. 8N 


1 5 5 . 9 W 


56 . 1 


34. 5 


28.8 


24. 1 


21.6 


5 4. eii 


155. 9W 


90. 1 


57. 8 


43.0 


26.6 


10.0 


3 0. 5N 


139 . 1W 


63.6 


56. 5 


50. 9 


41.1 


38.2 



The 


mean devi 


at ion (column 5 of 


Table II) 


generally 


increases with 


latitude up to 42. 7N 


but then 


decreases at 


the 


two highest 


latitudes. Additionally, the 


contribution 



of modes 2 to 5 in reducing the mean deviation is signifi- 
cant at the ncithern most points. To examine the relevance 
of these results, the cutoff period (minimum period) for the 
propagation cf free Rossfcy waves of mode n (see Emery and 
Magaard (1976)) given by T = (4irfX„)//j was evaluated. The 
results are shewn in Figure 26. The dotted line delineates 
the annual period and the dashed line the five year period. 
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Figure 26. Cutoff Period versus Latitude 
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The curves cf Figure 26 shew the following results: 



• Eased upon the theoretical cutoff oeriod annual oeriod 
taroclinic free Rossby waves may exist nor*h cf 25 n but 
only mods 1 can occur. 



• Earoclinic free Hossby waves at the annual period can- 
not exist north of approximately 36N (recall Magaard 
(1983) stated that first mode baroclinic Rossby waves 
were primarily evident between latitudes 30 to 40 North 
in the North pacific Ccean) . 

• Baroclinic free Rossby waves at the five year period 
cannet exist north cf approximately 52N, 

The results given above are significant. Since the free 
Rossby wave analysis was conducted using only ten years of 
model simulated temperature data, free Rossby waves of peri- 
ods greater than approximately five years can not be ade- 
quately resolved. Therefore, the results obtained for model 
grid point 54. 6N 155. 9W losses significance. Likewise, at 

grid point 48. 8N 155. 9W the theoretical cutoff period curves 
show that only firsn mode baroclinic Rossby waves an periods 
greater than approximately three years can exisn . There- 
fore, the low deviation at latitude 43. 3N resulning prima- 
rily from modes 2 to 5 is not necessarily due to baroclinic 
Rossby waves. The opposine may be said of the selecned grid 
points at latinudes 18.3N to 36. 6N. At these latitudes the 
first baroclinic mode is dominant and the cutoff period is 
4.5 months at 13. 3N and 12.8 months at latitude 36. 6N. At 
latitude 42. 7N the cutoff period is 21.4 months. However, 
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only the first and second modes can exist at periods less 
than five years. The mean deviation is high when only the 
first and second modes are fitted to the isotherm displace- 
ment fluctuations (96.2% and 88.1% in columns 1 and 2 , 
respectively, of Table II). Thus at model grid point 42. 7N 
155. 9W the model simulated isotherm displacemenr fluctua- 
tions cannot be explained by free baroclinic Rossby waves in 
the time domain. 

Thus far the examination of baroclinic free Rossby wave 
signature in the model simulated temperature field (or iso- 
therm displacement) has teen limited to the time domain 
only. In order to conduct an analysis of baroclinic free 
Rosshy wave signature in the space domain requires carrying 
out a cross-sprectral fit cf the x and y direction wave num- 
bers at the appropriate frequencies. In the context cf this 
study a cursory examination of the wave number vectors will 
be made. The locus of wave number vectors for mode 1 at 
appropriate frequencies for latitudes 18. 3N and 36. 6N are 
shown in Figures 27 and 28. The locus of wave numbers was 
determined using equation (11) for three different frequen- 
cies . 
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Figure 27. Locus of Wave Humber Vectors for Mode 1 ar 18. 3H 
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Figure 28. Locus of Wave Number Vectors for Mode 1 at 36. 6N 
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Ihs spacing between model grid points urilized in this 
study is cr tbs crder of 226 km in the north/south direction 



and 251 km (at 36N) in the east/west direction. The shor- 
test wavelength resolvable by the model is equal to twice 
the spacing between grid poinrs , which in this case is on 
the crder cf approximately 500 km. Therefore, Rossby waves 
resolvable by the model must have wavelengths greater than 
or equal to the 2ad value, where Ad is the grid spacing. 
The wave numbers in the x and y direction which correspond 
to a wavelength which is greater than the 2Ad wavelength are 
located withir the dotted rectangle in Figures 27 and 28. 
Stated another way, wave numbers outside of the boxed-in 
region correspond to wavelengths too small to be resolved by 
the mcdei. 



The inner circle in Figure 27 is the locus of wave num- 
bers found when a frequency corresponding to a six-month 
period is used. The middle circle and ourer crrcle are for 
frequencies ccrresponding to periods of nine months and 12 
months, respectively. Note that due to the different scale 
of eaca axis the circles appear as ellipses. Portions of 
all three circles fall within the enclosed region indicating 
that wavelengths greater than the 2nd wavelength can occur 
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at the af orement icned frequencies (periods). The circles in 
Figure 28 are icr frequencies corresponding to periods of 18 
months (inner circle), 24 months (middle circle) and 30 
months (cuter circle) , respectively. It is nored rhat at 
36. 6N only waves having a period of about 30 mcnrhs or 
greater have wavelengths long enough to be resolvable by the 
model. It follows that waves of 44d or 8ad would have even 
greater periods. Therefore, increasing latitude requires an 
increase in the Rcssby wave period to allow for wavelengths 
long enough to be resolved by the model grid spacing. Addi- 
tionally, using equation (11), it may be shown that for 
wavelengths that are resolved by the model grid, k^ + S^« An 
and thus the period 2n /a must be larger rhan the cutcff 
period. At high latitudes such long periods would not he 
well resolvad ty the ten-year simulation. Thus, even though 
low deviations result at high latirudes (column 5 of Table 
II) the isotherm displacemenr variability cannot be attrib- 
uted to actual Rossby waves with any degree of certainty. 
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IV. 



CONCLUSICNS AND RECOMMENDATIONS 



Eased upon the preceeding analysis the following results 



are summarized: 



• On the average, free baroclinic Rossby wave signature 
in the acdel simulated isotherm disolacement fluctua- 
tions is evident in the time domain* for the selected 
grid points south of 40N. 

• At the selected grid points south of 40N the first 
fcaroclinic mode shows predominance. 



• Eased upon the locus of wave number ellipses, Rossby 
waves at the appropriate periods are long enough to be 
resolvable witnic the model grid spacing. 

• Eased upon the deviation term^ certain times within the 
ten-veaf aodel run period indicate that the first five 
baroclinic modes do not adequately fit the isotherm 
displacement fluctuations (i.e., the deviation is high, 
therefore, the isotherm displacement fluctuations can- 
not be fullv explained by baroclinic free Rossby 
waves ) . 



• At selected grid points north of 40N the analysis 
results are inconclusive. 



• There is reasonably gccd agreement between the analysis 
results at model grid point 30. 5N 139. 1W and that con- 
ducted by Emery and Magaard (1976) for weather station 
November (3CN i40H) . 



The results stated in point number five above, are a product 



of several factors. The model run time (ten years) sets a 



constraint upon the periods that may be adequately resolved 



by the analysis method. The cutoff periods of the selected 



grid points ncrth of 40N gradually begin to approach a 



five-year period, and wavelengths resolved by the model grid 
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result in periods considerably greater than rhs curoff (min- 
imum) period. Periods greater than five years can nor be 
adequately resolved by the analysis method. A longer run 
time is required. Secondly, since the ocean model is a 
'•fresh'* water model the effect is to dimenish the magnitude 
of the Brunt — Vaisala frequency, particularly at higher 
latitudes. This results in eigenvalues that are large com- 
pared to "real'* ocean values { X is inversely proportional to 
N) . The larger eigenvalues then, in turn, increase the 
value of the cutoff period more than what occur when only 
the Coriolis parameter is considered. Finally, a seasonal 
signal, resulting from the seasonally varying model input 
parameters, dominates the isotherm displacement fluctuations 
over a greater portion of the water column at the northern 
selected grid points. Thus, the seasonal signal, which is 
primarily a thermodynamic effect, contaminates the first 
five modes. 



Eased upon these results several recommendations for 
further research may be set forth: 



• The seasonal signal resulting from thermodynamic 
effects should be removed. This could be achieved by 
either neglecting the model layers above the main ther- 
mccline or by applying a filter technique to remove the 
the regular seasonal signal. The former suggestion 
would oe impractical at the high latitude model grid 
points, however, due to the penetration of the seasonal 
signal down to the bottom of the water column. 
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• This study set the mean current equal to zero. Further 
studies should include the mean current in -^he analy- 
sis, particularly in the region of the California cur- 
rent, and North Pacific Current. A study by Kang and 
Maaaard (1979) has shown that the North Pacific Current 
has a significant impact on the first barcclinic shear 
mode . 



• Since the model neglects salinity further research 
should study the feasability of incorporating salinity. 
This would provide fcr more realistic density profiles 
and improve the analysis results at t ne higher 
latitudes . 
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• All further studies should include comparison to appro- 
priate observational studies when available, and 
thereby include a model evaluation. 

This study has teen a first step in analyzing the low fre- 
quency dynamics of the Haney (1980) ocean model below the 
mixed layer. A complete understanding of the dynamics of 
the large-scale variability of the ocean temperature field 
below the mixed layer (t hermocl ine) is the primary goal. 
This understanding will in turn lead to improvements in the 
ocean model and its diagnostic and prognostic capabilities. 



the ultimate goal, 
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